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ABSTRACT 

The presence of relativistic particles at the center of our galaxy is evidenced by the diffuse TeV 
emission detected from the inner ~2° of the Galaxy. Although it is not yet entirely clear whether 
the origin of the TeV photons is due to hadronic or leptonic interactions, the tight correlation of 
the intensity distribution with the distribution of molecular gas along the Galactic ridge strongly 
points to a pionic-decay process involving relativistic protons. In earlier work, we concluded that 
point-source candidates, such as the supermassive black hole Sagittarius A* (identified with the 
HESS source J1745-290), or the pulsar wind nebulae dispersed along the Galactic plane, could 
not account for the observed diffuse TeV emission from this region. Motivated by this result, we 
consider here the feasibility that the cosmic rays populating the Galactic Center (GC) region are 
accelerated in situ by magnetic turbulence. Our results indicate that even in a highly conductive 
environment, this mechanism is efficient enough to energize protons within the intercloud medium 
to the > TeV energies required to produce the HESS emission. 
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Observations of the Galactic Center (GC) with 
the High-Energy Stereoscopic System (HESS) 
have revealed the presence of diffuse TeV emis- 
sion spread out roughly 0.2° in Galactic latitude b 
within the inner 2° of our galaxy (Aharonian et al. 
2006). The strong correlation between the > 200 
GeV emission and the ~10* Mq of molecular gas 
distributed along the GC ridge, as traced by its 
CO and CS line emission (see, e.g., Tsuboi et al. 
1999), points to the decay of neutral pions pro- 
duced by the scattering of relativistic cosmic rays 
with the proton-rich target of overlapping clouds 
as the dominant source of this diffuse radiation 



(see, e.g., Crocker et al. 2005; Ballantyne et al. 
2007). 

The origin of these energetic hadrons is an in- 
triguing puzzle because the observed gamma-ray 
spectrum requires an underlying cosmic-ray pop- 
ulation quite different from that seen at Earth. 
Specifically, the gamma-ray spectrum measured 
by HESS in the region |/| < 0.8° and \b\ < 0.3° 
(with point-source emission subtracted) can be 
reasonably fit with a power law with photon in- 
dex r = 2.29 ± 0.27. Since the spectral index of 
the gamma rays tracks the spectral index of the 
cosmic rays themselves, the implied cosmic ray in- 
dex (~2.3) is then much harder than that (~2.75) 
measured locally. 
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Possible sources of energetie hadrons at the GC 
were recently considered by Wommer et al. (2008). 
The results of this effort seemingly rule out point 
sources such gittarius A* and the pulsar wind 
nebulae dispersed along the Galactic ridge, and 
thereby give credence to the possibility that the 
relativistic protons arc accelerated throughout the 
GC medium. Following up on this result, Fatuzzo 
& Melia (2011) found that stochastic accelera- 
tion by magnetic turbulence within the intcrcloud 
medium, which is effectively a 1-D random walk 
in energy process, can produce a distribution of 
particle energies whose high-energy tail is capable 
of reproducing the HESS data, so long as the tail 
extends to > 1 TeV energies. 

However, the electric field used to calculate the 
energy evolution of particles in Fatuzzo & Melia 
(2011), which was derived directly from the form 
of the turbulent magnetic field using Faraday's 
Law, had a component parallel to the overall mag- 
netic field. As a result, the electric field was so effi- 
cient at energizing protons that the rciquircd parti- 
cle distributions could only be produced if the tur- 
bulent field was much weaker than the underlying 
field or if acceleration was limited to small "active 
regions" within the intercloud medium. However, 
the high conductivity of the medium within the 
GC environment makes it highly likely that a com- 
ponent of the electric field parallel to the underly- 
ing magnetic field would be quickly "quenched" . 

The purpose of this paper is to reassess the fea- 
sibility of stochastic acceleration within the GC 
region with the (more realistic) assumption that 
the highly conductive medium does not allow for 
an electric field component parallel to the mag- 
netic field. Toward that end, we adopt the for- 
malism of O'SuUivan ct al. (2009) to construct 
the turbulent fields so that the electric field is ev- 
erywhere perpendicular to the total magnetic field. 
The spatial and energy diffusion of cosmic-ray pro- 
tons within the molecular cloud and intercloud re- 
gions are then investigated via numerical simula- 
tions. Specifically, both the spatial and energy 
diffusion coefficients over a relevant range of pa- 
rameter space are calculated and used to compare 
estimates of the time required to energize protons 
up to TeV energies with the escape and cooling 
times associated with both the cloud and inter- 
cloud environments. 

Our results indicate that protons in the inter- 



cloud medium can be cmcrgizcd up to the > TcV 
energies required to produce the observed HESS 
emission. As such, stochastic particle accelera- 
tion by magnetic turbulence appears to be a viable 
mechanism for Cosmic-ray production at the GC. 

2. The Physical Conditions 

The large concentration (up to ^10^ -^0) of 
dense molecular gas at the GC is largely con- 
fined to CMC's with a size ~50-70 pc (Glisten & 
Philipp 2004). These clouds appciar to be clumpy 
with high-density (^10^ cm^'^) regions embedded 
within less dense (~10'^ '' cm~^) envelopes (e.g., 
Walmsley et al. 1986) and are threaded by a 
pervasive magnetic field whose milligauss strength 
is suggested by the rigidity of non-thermal fila- 
ments interacting with the molecular clouds (see, 
e.g., Yusef-Zadeh & Morris 1987; Morris & Yusef- 
Zadeh 1989; Morris 2007). Confirming evidence 
for such field strengths in and around the CMC's 
is provided by their apparent stability. The ob- 
served pressure Ppiasma~ 10~^'^ erg cm~^ due to 
the hot plasma between the clouds is an order 
of magnitude smaller than the turbulent pressure 
P(„rf)~10^* erg cm~^ within the GMC environ- 
ment (Glisten & Philipp 2004), seemingly ruling 
out pressure confinement. If clouds are instead 
bound by their own magnetic fields, then equat- 
ing the turbulent and mag netic (SVStt) energy 
densities gives field strengths of ^0.5 mC within 
the clouds, not too different from the typical value 
measured in the non-thermal filaments. 

We are now also reasonably sure of the mag- 
netic field strength between the clouds. In the 
past, the field intensity near the GC had been un- 
certain by two orders of magnitude. We've just 
seen how on a scale of ~100 pc, field strengths 
can be as high as -^1 mG. At the other extreme, 
equipartition arguments based on radio observa- 
tions favor fields of only ^6 on ^400 pc 
scales (LaRosa et al. 2005). But a more care- 
ful analysis of the diffuse emission from the cen- 
tral bulge has revealed a down-break in its non- 
thermal radio spectrum, attributable to a tran- 
sition from bremsstrahlung to synchrotron cool- 
ing of the in situ cosmic-ray electron population. 
Crocker et al. (2010) have shown recently that 
this spectral break requires a field of ^^50 /iG 
extending over several hundred parsecs, lest the 
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synchrotron-emitting electrons produce too much 
7-ray emission given existing constraints (Hunter 
et al. 1997). 

While the structure of this magnetic field is 
not well understood, magnetic fluctuations are ex- 
pected to be present in essentially all regions of 
the interstellar medium. For example, molecu- 
lar clouds are observed to have substantial non- 
thermal contributions to the observed molecular 
line-widths (e.g., Larson 1981; Myers, Ladd, & 
Fuller 1991; Myers & Gammie 1999). These non- 
thermal motions are generally interpreted as aris- 
ing from MHD turbulence (e.g., Arons & Max 
1975; Gammie & Ostriker 1996; for further evi- 
dence that the observed linewidths are magnetic 
in origin, see Mouschovias & Psaltis 1995). In- 
deed, the size of these non-thermal motions, as in- 
dicated by the observed line- widths, are consistent 
with the magnitude of the Alfvcn speed (e.g., My- 
ers & Goodman 1988; Crutcher 1998, 1999; Mc- 
Kee & Zweibel 1995; Fatuzzo & Adams 1993). As 
a result, the fluctuations are often comparable in 
magnitude to the mean values of the fields. 

For this work, we treat molecular clouds within 
the GC as spherical (i?c = 30 pc), uniform density 
(riH = n^f2/2 = 10'* cm^'^) structures threaded by 
an underlying uniform magnetic field Bq = BqZ, 
where Bq = 0.50 mG. The Alfven speed within 
the cloud environment is therefore taken to be 
PS 11 km/s. Consistent with the limits placed 
by Crocker et al. (2010), we treat the intercloud 
medium spherical (Ric = 200 pc) , low density 
{uh = 10 cm~^) structure threaded by an un- 
derlying uniform magnetic field Bq = Bqz, where 
Bq = 50 /iG. The Alfvcn speed within the inter- 
cloud environment is therefore va = 35 km/s. For 
both regions, we assume that the magnetic turbu- 
lence has the same energy density as that of the 
underlying uniform magnetic field. 

For completeness, we note that the conditions 
much closer to Sagittarius A* are somewhat differ- 
ent and appear to be controlled primarily by ongo- 
ing stellar wind activity (Rockefeller et al. 2004). 
But this is a very small region compared to the 
rest of the TeV emitting gas, so we do not expect 
it to significantly infiuence our results. 



3. The Turbulent Fields 

The standard numerical approach for analyz- 
ing the fundamental physics of ionic motion in 
a turbulent magnetic field treats the total mag- 
netic field B as a spatially fluctuating component 
SB superimposed onto a static background compo- 
nent Bq, where SB is generated by summing over 
a large number of randomly polarized transverse 
waves with wavelengths A„ = 2n/kn, logarithmi- 
cally spaced between Xmin and Xmax (e.g., Gian- 
coh & Jokipii 1994; Casse et al. 2006; O'Sullivan 
et al. 2009; Fatuzzo et al. 2010). Adopting a 
static turbulent field removes the necessity of spec- 
ifying a dispersion relation between the wavevec- 
tors kn and their corresponding angular frequen- 
cies LJn- This approach therefore allows one to con- 
sider highly non-linear turbulence {SB >> Bq), 
or even remove the background component alto- 
gether. Of course, turbulent magnetic fields in 
cosmic environments are not static. Nevertheless, 
a static formalism in spatial diffusion calculations 
of relativistic particles is justified for environments 
in which the Alfven speed is much smaller than the 
speed of light. 

This paper focuses on the energy diffusion of 
cosmic rays propagating through a turbulent mag- 
netic environment, which then requires the use 
of a time-dependent formalism in order to self- 
consistently include the fluctuating electric fields 
that must also be present. Toward that end, we 
assume that the GC environment is well repre- 
sented by a nonviscous, perfectly conducting fluid 
threaded by a uniform static field Bq = Bqz, 
and use linear MHD theory to guide us. In gen- 
eral, three types of MHD waves exist in the lin- 
ear regime — Alfven, fast and slow. Following the 
formalism of O'Sullivan et al. (2009), we consider 
here only Alfven waves, so that the turbulent mag- 
netic field is defined by the sum of N randomly 
directed waves 

N 

SB = J2^n e'(^"-^"-'^-*+^") , (1) 

n=l 

where the direction of each propagation vector fc„ 
is set through a random choice of polar angles 6n 
and (pn , and the phase of each term is set through 
a random choice of 

Alfven waves don't compress the fiuid through 
which they propagate, and are therefore charac- 
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terized by a fluid velocity v that satisfies the; con- 
dition k ■ V = 0. This condition in turn implies 
that V ■ Bo = 0. We can therefore write the fluid 
velocity associated with the n-th term in Equation 
(1) as 

Svr. = ±A„ ^ 1^ ,iik.-r-^nt+p.) ^ (2) 
■Do \z X k„\ 

where the sign is chosen randomly for each term 
in the sum. The dispersion relation for Alfvehic 
waves is given by the expression 



w„ = ?;^fc„| cos6'„| , 



(3) 



where va is the Alfven speed, and 0„ is the angle 
between fc„ and Bq 

The corresponding magnetic field for each wave 
then follows from the linear form of Ampere's Law, 
as given by 



' z z x kj^ 



(4) 



This formalism is identical to that of O' Sullivan et 
al. (2009), with the exception that we allow for a 
random choice of sign in each term. We find, how- 
ever, that including a random sign has no effect 
on the statistical output measures of our numeri- 
cal simulations. 

The total number of terms in the sum is given 
max/kmin]^ and tlic vaiucs of k^ 
arc evenly spaced on a logarithmic scale between 

kmin = 27r/A„,a3; and k^ax ^ '^T^/^min- TllC dc- 

sired spectrum of the turbulent magnetic field is 
set through the appropriate choice of F in the scal- 
ing 



— ^1 



1 -r 



Afc„ 



-r+i 



(5) 



(e.g., r = 3/2 for Kraichnan and 5/3 for Kol- 
mogorov turbulence). The value of Ai is set by 
a parameter ^ that specifies the averaged energy 
density of the turbulent field via the definition 



{SB') = {Y^^dB^.5B:,) 



n n' 

-r+1 



kn 



(6) 



where the n ^ n' terms average to zero. We note 
that for our adopted scheme, the value of /S.kn/kn 
is the same for all values of n. We further note that 
^ = 2 corresponds to the real part of the turbulent 
field having the same energy density as a uniform 
field So. 

Naively extending the results from linear MHD 
theory to our formalism, one would obtain the to- 
tal electric field 5E associated with the turbulent 
magnetic field defined by Eq. (1) by summing over 
the terms 

5 En =-6VnX Bo. (7) 

Although SE - Bo — 0, the second order term 
SE-SB 0. The presence of an electric field com- 
ponent parallel to the magnetic field in this second 
order term can significantly increase the accelera- 
tion efficiency artificially, especially if the formal- 
ism is extended to the nonlinear regime {SB ^Bo). 
However, the interstellar medium is highly con- 
ductive, and as such, any electric field component 
parallel to the magnetic field should be quickly 
quenched. To circumvent this problem, we adopt 
the formalism of O 'Sullivan et al. (2009) and first 
obtain the total fluid velocity 5v via the summa- 
tion 



N 



Sv = SVn 



(8) 



We then use the MHD condition to set the total 
electric held: 



SE = J-^xB, 
c 



(9) 



where B = Bqz + SB. 

4. Numerical Analysis 

The equations that govern the motion of a 
relativistic charged proton with Lorentz factor 7 
through the turbulent medium are 



and 



df 
'di 



V X B 
SE-\ 



= V 



(10) 



(11) 



Although these equations are deterministic, the 
chaotic nature of motion through turbulent flelds 
necessitates a statistical analysis. 
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We define a single experiment as a numerical 
investigation of particle dynamics through a given 
environment and a given particle injection energy. 
The environment is specified by the underlying 
field strength Bq and Alfven speed va (see §2), and 
the turbulent fields are specified by the parameters 
r, ^rnax, ^mim Nk and ^. For cach experiment, we 
numerically integrate the equations of motion for 
Np = 1000 protons randomly injected from the 
origin with the same initial energy, as specified by 
the Lorentz factor 70. The equations of motion 
are integrated for a time At = lOOXmax/c, with 
each particle sampling its own unique magnetic 
field structure (i.e., the values of /3„, and 
the choice of a ± are chosen randomly for each 
particle). 

As is well known, the diffusion coefficients pro- 
vide a useful output measure for the characteri- 
zation of the diffusion process since their values 
are constant once the particles are in the diffusion 
regime. We therefore adopt = (A7^)/(2Ai), 
= (A,T2)/(2At), and £>|| = {Az'^)/{2At) as 
the output measures of our experiments (recall 
that the underlying field Bq is in the z direction). 
The diffusion constants can then be used to ob- 
tain the "acceleration time" Tacc = lo/^i^ which 
characterizes how long it would take low energy 
particles to diffuse to •jonipC^ energies (and hence, 
attain a high energy tail > 7oTOpC^), and the es- 
cape time Tesc = /D\\, which characterizes how 
long it would take those particles to diffuse a dis- 
tance R along the preferential z direction. 

As elaborated upon in §3 below, a given ex- 
periment is defined by the parameters Bq, va, T, 
Xmaxi s-nd 70, and the output measures are the 
diffusion coefficients D^, D± and £)||. These val- 
ues are summarized in Table 1 for all experiments 
performed in this work. 

4.1. Baseline experiment 

Our first goal is to find a numerical scheme that 
minimizes computer time without sacrificing accu- 
racy. Toward that end, we perform our first nu- 
merical experiment for 70 = 10^ particles injected 
into the intercloud environment {Bq = 50 /iG, 
Va = 35 km/s). We adopt a baseline set of tur- 
bulent field parameters F = 5/3, X„iax = 1 pc, 
A„i„ = 0.002 pc, C = 2 and Nk = 25. We note 



that the particle radius of gyration, 

^.o = 0.02pc(^)(^)"\ (12) 
falls comfortably within the values of Xmin and 

Xmax ■ 

As can be seen from Figure 1, a single particle's 
energy, as characterized by A7/70 = 7/70 ^ Ij 
changes in a random-like fashion. As such, the 
energy distribution for an ensemble of particles 
injected with the same energy becomes normal 
once the particles have fully sampled the turbu- 
lent nature of the accelerating electric fields. This 
point is clearly ilhistrated by Figiire 2, which 
shows the distribution of A7 values at time t = 
Wi)X,nax/c for the 1000 particles tracked in our 
baseline experiment (Experiment 1). One can 
therefore quantify the stochastic acceleration of 
particles in turbulent fields through the variance 
of the resulting distributions of initially mono- 
energetic particles. To illustrate this point, we plot 
in Figure 3 the variance cr^ = ^/{A^ of the A7 
distribution function of time for the particles 
tracked in our baseline experiment. As expected 
from the random nature of stochastic acceleration, 
a-y oc \/i once particles have had a change to sam- 
ple the turbulent nature of the underlying fields, 

i.e., for t > Xmax I c. 

The spatial diffusion of particles will be differ- 
ent in the parallel and perpendicular directions to 
the underlying magnetic field Bqz^ resulting in dif- 
ferent variances in the distributions of particle dis- 
placement along and across the underlying field. 
To illustrate this point, we plot Gx = 1/ (Aa;^) and 
= ^ (Az^) as a function of time in Figure 4 
for the particles in our baseline experiment. As 
expected, particles diffuse farther in the direction 
parallel to the magnetic field, and both output 
measures become oc at time t > Xmax I c- 

A fundamental issue in this analysis is what 
value of TVfc will allow our discrete treatment of 
the turbulent field to adequately represent the 
continuous fields found in nature. Toward that 
end, we repeated our baseline experiment with 
Nk = 250. The difference in output measures 
were smaller than 10%, indicating that setting 
iVfc = 25 provides good accuracy in our results. 
We also repeated our baseline experiment with 
Xmin = 0.0002 pc. Again, the difference in output 
measures were smaller than 10%, indicating that 
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Fig. 1. — The fractional change in particle energy 
A7/70 as a function of time for a 70 = 10^ par- 
ticle injected into an intercloud-like environment 
{Bq = 50 /xG, VA = 35 km/s) with an Alfvenic tur- 
bulent field defined by the parameters \max — 1 
PC A„ji„ = 0.002 pc, r = 3/2, ^ = 2.0, and 
Nk = 25. 
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Fig. 3. — The variance of the particle energy 

distribution as a function of time for the Np = 
1,000 particles in our baseline numerical experi- 
ment (Experiment 1). The dotted line serves as a 
reference and has a slope of 1/2, clearly indicating 
that cr^ oc for time t > Xmax/c 



150 



100 




Fig. 2. — The particle energy distribution for 

the ensemble of Np = 1,000 particles at time 
t = lOOXmax/c in our baseline numerical exper- 
iment (Experiment 1). The solid line shows a 
Gaussian fit to the data. 
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Fig. 4. — The variance of the particle spatial dis- 
tributions perpendicular (x) and parallel (z) to the 
underlying field Bq as a function of time for the 
Np = 1, 000 particles in our baseline numerical ex- 
periment (Experiment 1). The dotted line serves 
reference and has a slope of 1/2, clearly indi- 
cating that (jj, and cr^ are proportional to for 
time t > Xmax/c. 
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our output measures are not sensitive to the value 
of Xmin so long as the particle radius of gyration 
Rg > Xmin (see also Fatuzzo et al. 2010). 

Guided by the results of our baseline analysis, 
we will perform the remainder of our experiments 
using Nk = 25 and Xmin = O.lRgo, where 



R. 



eBo 



(13) 



is the radius of gyration for the injected particles 
in the absence of turbulence. The; governing equa- 
tions for each particle will be integrated out to a 
time sufficiently long to ensure that and (7z 

arc proportional to \/t by the end of the integra- 
tion (usually At = lOOXmax/c), and the particles 
have therefore fully sampled the turbulent nature 
of the magnetic and electric fields. This scheme is 
expected to provide accurate results with minimal 
computing resources. 

As a consistency check, we perform a set of ex- 
periments (2 - 7) using the same physical condi- 
tions corresponding to models W3 and X3 pre- 
sented in O'Sullivan et al. (2009; see Table 1 
and Figures 3 and 4). Specifically, we consider 
a physical environment defined by the parameters 
Bo = 3 fxG and va = 0.002 c, and a turbulent field 
defined by F = 5/3 (Kolmogorov) and Xmax = 1 
kpc for ^ = 0.2 (corresponding to model W3 for 
which {6B/Bo)^ = 0.1) and ^ = 2 (corresponding 
to model X3 for which {SB/Bq}'^ = 1). We com- 
pare our results for the acceleration time Tacc 

to 

those obtained by O'Sullivan et al. (2009) in Fig- 
ure 5. While there is general agreement between 
these results, our values of the acceleration times 
are consistently about a factor of two greater than 
those to which we are comparing. 

As a final consistency check, we compare in 
Figure 6 the spatial diffusion coefficients D± and 
D|| obtained for Exp. 5 7 to those obtained 
by Fatuzzo et al. (2010) under the same physical 
conditions. We note that the turbulent magnetic 
field used in our earlier work is of the form given 
by Equation 1, but with 



An = An (cosa„ y' ± i sina„ z') , 



(14) 



where the y' — z' plane; is normal to the k direction 
and an is picked at random for each n term. As 
can be seen from Figure 6, there is good agree- 
ment between results, indicating that the turbu- 



lent magnetic field adopted in this work is equiva- 
lent to that used in our earlier work. We note also 
that the presence of a weak turbulent electric field 
(i.e., SE << SB), which was not included in the 
calculations of Fatuzzo et al. (2010), has a negli- 
gible effect on the spatial diffusion of particles. 

4.2. Survey of Parameter Space 

As noted above, the spatial and energy diffu- 
sion of particles through turbulent fields is not 
sensitive to the minimum wavelength so long as 
the particle's radius of gyration exceeds X^in- For 
a given environment, as defined by Bo and va, the 
diffusion process is thus dependent upon the max- 
imum turbulence wavelength Xmax, the turbulent 
field strength, as characterized by ^, and the tur- 
bulence spectrum, as characterized by the spectral 
index F. 

Quasi- linear theory predicts that the energy dif- 
fusion coefficient for relativistic particles scales as 



^ ~ c2 B^o 



Rn 



Xn 



r-1 2 2 

RgC ' 



(15) 



where p is the particle momentum (Schlickeiser 
1989). For relativistic particles, the energy diffu- 
sion coefficient should therefore scale to our model 
parameters as 



a ^ A, 



max '0 ' 



(16) 



Equation (15) however does not appear to be valid 
in the strong turbulence limit (O'Sullivan ct al. 
2009). We therefore investigate how the energy 
diffusion coefficient depends upon Xmax and ^ for 
Kolmogorov (F = 5/3) turbulence. Specifically, 
we perform two sets of experiments designed to 
explore parameter space around our baseline val- 
ues. For the first set (Experiments 8 - 10), we 
vary Xmax from 0.32 - 10 pc, while in the second 
set (Experiments 11 - 13) we vary ^ between 0.2 
^ 6.4. 

The energy diffusion coefficients for Experi- 
ments 1 and 8-10 are presented in Figure 7, 
with the best line fit to the data indicating that 



A; 



-0.47 



This results confirms that quasi- 



linear theory, which predicts that cx A„ 



for 



Kolmogorov turbulence, is not applicable in the 
strong turbulence (^ > 1) limit. The energy diffu- 
sion coefficients for Experiments 1 and 8-10 are 
presented in Figure 8, with the best line fit to the 
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Fig. 5. — Acceleration time Tacc as a function of 
the injected particle Lorentz factor, for both our 
work and the work of O'SuUivan ct al. (2009). 
Open squares (W3) and circles (X3) denote the 
results of O'SuUivan et al. (2009). Solid squares 
and circles denote our results. We note that Ac = 
0.77A„„^/(27r). 
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Fig. 6. — Diffusion coefficients D±_ and D\\ as a 
function of rigidity p = 2TTRgQ/\„iax for experi- 
ments 5-7 [solid squares and circles] compared 
with results obtained by Fatuzzo et al. (2010) for 
which ^ = 2 (open squares and circles). 



data indicating that D-^ oc ^^-^ in both the weak 
and strong turbulence limit. 

5. Application to the GC Environment 

The stochastic acceleration of protons in a mag- 
netically turbulent environment is essentially a 1- 
D random walk process that results in a particle 
energy distribution that broadens in time. Our 
goal in this section is to calculate the characteris- 
tic particle energy for both the inter-cloud medium 
and the molecular cloud region observed at the 
GC. As shown in §4, a numerical treatment of 
particle diffusion requires an integration time-step 
St ~0.1Ami„ ^O.OlRg, but a total integration time 
that is At ~100Amox- Numerically integrating the 
equations of motion for protons that are energized 
from a thermal state to TcV energies is therefore 
not computationally feasible given the small ra- 
dius of gyration of thermal particles. 

As such, we determine the characteristic parti- 
cle energy in each region by comparing the "accel- 
eration time" Tacc = lo/Dj, which characterizes 
how long it would take a distribution of low-energy 
particles to attain energies ^jQinpC^, with the es- 
cape time Tesc = which characterizes how 
long it would take those particles to diffuse a dis- 
tance R along the preferential z direction, and the 
proton cooling time Tpp, which characterizes how 
long a proton can move through a region before 
losing a significant fraction of its energy to pp scat- 
tering with the ambient medium. 

To determine how the acceleration and escape 
times depend on particle energy for the intcr-cloud 
medium {Bq = 50 /zG, va = 35 km/s, R = 200 
pc), we perform a set of Experiments (14 - 18) 
using our base-line parameters T = 5/3, Xmax = 1 
pc, and ^ = 2. We perform two additional experi- 
ments (19 & 20) at 70 = 10^ and 70 = 10*^ for the 
same environment, but with F ~ 3/2. We perform 
a final set of Experiments (21 - 26) to determine 
how the acceleration and escape times depend on 
particle energy for the molecular cloud medium 
{Bo = 500 iiG, VA = 11 km/s, i? = 30 pc). 

The energy loss rate of protons with energies 
needed to produce 7r°-decay 7-rays is dominated 
by nuclear energy losses due to pp scattering with 
the ambient medium (Aharonian & Atoyan 1996). 
As such, the cooling time, Tpp, of the protons de- 
pends on the pp-scattering cross-section, app, and 
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Fig. 7. — The energy diffusion coefficients Dj as a 

function of maximum turbulence wavelength A„,„ 
for experiments 1 and 8 - 10. The best line fit to 
the data has a slope of —0.47. 
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Fig. 8. — The energy diffusion coefficients as a 
function of ^ for experiments 1 and 11 — 14. The 
best line fit to the data has a slope of 1.2. 



the inelasticity parameter, k. Over a broad range 
of proton energies, neither of these quantities sig- 
nificantly varies so the usual method is to adopt 
the constant average values cr-pp « 40 mb and k ss 
0.45 (see, e.g., Markoff et al. 1997). That being 
the case, the proton cooling time becomes inde- 
pendent of proton energy: 

Tpp = {nHCKapp)~^ , (17) 

where Uh is the number density of ambient pro- 
tons. The value of Tpp is therefore s::^ 2 x 10^'* s for 
the inter-cloud region and « 2 x 10^^ within the 
molecular cloud region. 

Our principal results are presented in Figures 

9 and 10. A characteristic particle energy for 
each region can be estimated by the crossing point 
where Tqcc becomes greater than either Tesc or Tpp. 
Interestingly, this value is ^3 Tev for the inter- 
cloud region, indicating that the TeV cosmic-rays 
associated with the HESS observations can thus 
be produced via stochastic acceleration within the 
turbulent intercloud environment. This result ap- 
pears to hold true for both Kolmogorov (F = 5/3) 
and Kraichnan (F = 3/2) turbulence, and suggests 
that energy diffusion in the strong turbulence limit 
is not sensitive to the turbulence power-spectrum. 
In contrast, while the acceleration time associated 
with the molecular clouds are comparable to those 
in the intercloud medium, the escape time and the 
proton cooling time are significantly shorter. As 
such, it is clear from Figure 10 that TeV protons 
cannot be energized within the molecular cloud 
environment. 

6. Conclusion 

The purpose of this paper has been to assess the 
feasibility of stochastic acceleration within the GC 
region in producing the TeV cosmic-rays revealed 
by the diffuse HESS emission correlated with the 
molecular gas distributed along the GC ridge. As 
shown in previous work, this emission can be pro- 
duced by cosmic rays with an energy distribution 
that has a high-energy tail extending out beyond 
a few TeV. Such a distribution is naturally pro- 
duced by the stochastic acceleration of particles in 
a magnetically turbulent environment (which is ef- 
fectively a 1-D random walk in energy). Thus, the 
remaining question is whether or not this mecha- 
nism is efficient enough to energize protons to the 
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Fig. 9. — The acceleration time Tacc and escape 
time Tesc for the intercloud environment as a func- 
tion of injection energy (as characterized by the 
Lorcntz factor 70) for experiments 14 20 (solid 
squares and circles) and experiments 19 - 20 (open 
squares). The solid lines represent the best-line 
fits to the data from experiments 14 20, and the 
dashed line represents the value of Tpp in the in- 
tercloud region . 
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Fig. 10. — The acceleration time Tacc and escape 
time Tesc for the molecular cloud environment as 
a function of injection energy (as characterized by 
the Lorentz factor 70) for experiments 21 - 26. 
The solid lines represent the best-line fits to the 
data, and the dashed line represents the value of 
Tot, in the GC molecular clouds. 



required TciV cinergicis in the highly conductive in- 
terstellar medium permeating the galactic center. 

To resolve this issue, we performed a series of 
numerical experiments in order to calculate the 
spatial and energy diffusion coefficients for protons 
in both the molecular cloud and the intercloud 
medium located within the inner several hundred 
parsecs of the galaxy. We then estimated the char- 
acteristic energy of the resulting particle distribu- 
tion by comparing the acceleration time required 
for particles to reach a certain energy, the escape 
time which characterizes how long it takes protons 
to diffuse out of each region, and the cooling time 
which characterizes how long a proton can move 
through a region before losing a significant frac- 
tion of its energy to pp scattering with the ambient 
medium. 

Our results indicate that for the physical condi- 
tions observed in the intercloud medium, together 
with reasonable estimates of the turbulent field 
strength and maximum turbulence wavelength, 
protons can be accelerated to characteristic en- 
ergies ^3 TeV, indicating that the TeV cosmic- 
rays observed by HESS can thus be produced via 
stochastic acceleration within the turbulent inter- 
cloud environment. In contrast, while the accel- 
eration time associated with the molecular clouds 
are comparable to those in the intercloud medium, 
the escape time and the proton cooling time are 
significantly shorter. As such, protons cannot be 
energized to TeV energies within the molecular 
cloud environment. 

These results are very encouraging, not only be- 
cause the idea of stochastic acceleration in a tur- 
bulent magnetic field at the GC appears to be a 
viable mechanism for producing the cosmic rays 
observed in that region, but especially because the 
characteristic energy attained by the relativistic 
protons matches the observations very well. We 
note, in this regard, that the physical conditions 
used in our simulations are uniqiic to the GC. As 
such, cosmic rays like those observed by HESS are 
not easily produced by this mechanism anywhere 
else in the Galaxy. 

So the investigation now turns to the very im- 
portant subsequent question, which we already al- 
luded to in the introduction, viz. can we under- 
stand the origin of cosmic rays observed at Earth, 
at least up to energies -^3-5 TeV, as the result 
of stochastic acceleration at the GC followed by 
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energy-dependent diffusion and escape across tlie 
Milky Way? Simulations designed to address this 
issue are currently underway, and the results will 
be reported elsewhere. 
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Table 1 
Summary of Experiments 
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